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INTRODUCTION 


At  the  Air  Force  Geophysics 
Laboratory! AFCL)  Weather  Test  Fac i 1 i ty( WTF) 
at  Otis  Air  Force  base,  MA,  an  instrumented 
tower  network  has  been  utilized  to  explore 
techniques  for  the  measurement  of  slant 
visual  range  (SVR)  and  the  short  range  pre¬ 
diction  of  below-llmit  SVR  conditions.  SVR  is 
defined  as  the  slant  distance  to  the  farthest 
high  intensity  runway  edge  light  or  approach 
light  which  a  pilot  can  see  at  decision  height 
(DM)  on  the  approach  path.  This  paper  will 
present  the  results  ot  the  AFCL  development 
an  i  t  valuation  of  a  remote  tower  SVR  system 
under  Categories  I  and  II  operations.  Three 
prediction  techniques  (Markov,  KKKP  and 
Equivalent  Markov)  to  yield  probability 
estimates  of  be  low- limit  SVR  conditions 
are  evaluated  and  comparisons  made  to  deter¬ 
mine  their  respective  accuracy  and  reliability. 
Forecast  times  examined  are  2,5,10,30  and  bO 
mi  nutes • 

2.  TEST  FACILITY  AND  DATA  SETS 

Development  and  testing  of  a  SVR 
syst  n  are  a  part  of  the  continuous  evaluation 
ot  various  meteorological  measurement  in¬ 
struments  at  the  AFGL  WTF.  Measurements  ol 
atmospheric  extinction  coefficient  are  made 
by  an  array  of  E(»AG  forward  scatter  meters  (FSM) 
mounted  on  three  towers  (A,Q,  and  X)  in  the  WTF, 
f  fgure  1 • 

The  WTF  Is  located  in  the  Cape 
Cod  area  where  low  visibility  episodes  arc 
predominantly  caused  by  heavy  advecti on  fog 
accompanied  by  light  rain  or  drizzle.  The 
continuous  data  stream  from  August  1977  to 
April  1979  revealed  a  systematic  difference 
between  the  SVR  and  runway  visual  range  (RVR) 
measurements*  Figure  2  depicts  a  pro-1  ronl.il 
band  of  showers  tol lowed  by  a  period  of 
heavy  advectlon  fog.  Note  that  during 
periods  oi  moderate  rain,  the  vertical  gradient 
of  extinction  coefficient  (lO~lkm~I)  is 
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nearly  zero.  During  the  advection  peri, 
however,  there  is  a  significant  increase 
extinction  with  height. 


Figure  i.  Conf igural ion  of  Instrument  1 >  a\ 


Figure  2.  Tower  A  Extinction  Coefficients 

14-15  October  1978 

SVR  specification  equ.it  ions  de¬ 
veloped  from  a  single  set  of  these  two 
of  low  visibility  episodes  could  rediu  e  t  l,i* 
probability  of  detecting  a  be  low-limit  .  .  nd  t 
of  SVR .  In  order  to  remove  this  svMi* it  i « 
bias,  the  selection  of  low  vi si  hi  1  it v 
episodes  was  restricted  to  the  advent  ion  tog 
type.  Episodes  were  selected  bv  applying 
the  criterion  that  the  SVR  in  the  landing  .*«•! 
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(toner  A)  was  1km"1  or  more  for  *  period  of  on 
hour  or  greater*  Sped  Mention  end  prediction 
equations  were  defined  from  a  dependent  data 
set  which  consisted  of  twenty  reduced 
visibility  episodes  and  contained  in  excess  of 
6300  minutes  of  da:a.  The  independent  data  set, 
used  to  test  and  compare  the  various  techniques! 
was  drawn  from  twenty-five  other  episodes 
totalling  over  9900  minutes  of  reduced 
visibility. 

3*  SVR  SPECIFICATION 

3.1  Category  II 

Haring  and  Gelsler  (1978)  de¬ 
monstrated  the  accuracy  of  FSM  measurements 
as  specifiers  of  SVR  under  Category  11  conditions. 
Those  results  also  Indicated  that  during 
advectlon  fog  conditions  observed  at  Cape  Cod, 
a  30  ft  remote  tower  system  yields  SVR  estimates 
nearly  as  accurate  as  a  100  ft  remote  tower  system 
In  which  point  visibility  measurements  are  made 
at  the  same  height  as  the  Category  II  OH. 

The  earlier  study  compared  a 
method  which  used  remote  tower  measurements 
of  visibility  to  specify  approach  sone  SVR 
(ASVR100)  to  a  method  which  relies  on  the  touch¬ 
down  RVR  measurement.  Table  1  summarises 
the  relevant  statistics  for  the  tower  option 
which  relies  on  measurements  at  50  ft  CASVR 1 00 ) 
and  the  control  technique. 


STATISTIC 

/tfviSb 

RVR 

CORRCLATION  came  »f  NT 

.974 

.tu 

21.4 

4a  tf 

HAS 

14 

37  9 

ThR!  SHOUF  Uml 

TS 

to 

71 

POD 

ts 

71 

FAS 

3 

4 

THSCSHOtP  12.  im*1 

TS 

12 

47 

POD 

♦2 

47 

FAR 

0 

Table  1:  STATISTICS  CATEGORY  II  SPECIFICATION 

The  threshold  5  km"l  corresponds 
to  about  1/2  mile  (800  m)  daytime  visibility 
and  1  mile  (1600  m)  at  night.  The  threshold 
12  km*'  corresponds  to  1/6  mile  (600  m) 
daytime  and  1/2  mile  at  night*  Clearly, 

X5VRI03  yields  superior  spec  1 f icat ions  of 
ASVR100  chan  does  RVR.  Figure  3,  which  is 
a  time  scries  plot  of  an  advectlon  fog  episode, 
demonstrates  that  the  major  deficiency  of  the 
RVP,  method  lies  In  the  persistent  optimism  it 
conveys  to  the  pilot  st  DH  of  seeing  his  re¬ 
ference  point* 


3.2  Category  X 

Specification  equations  for  200 
ft  SVR  (ASVR200)  were  developed  using  FSM  measure¬ 
ments  st  discrete  points  In  the  vertical  are 
converted  to  a  weighted  vertical  average  which 
is  used  to  represent  SVR  in  the  approach  zone 
through 


Application  ol  a  multiple  *i"t*ar 
regression  technique  to  a  dependent  data 
yielded  specification  equations  for  A SV92.V  a* 
foil ows i 

EST1  *  Q100  ♦  .86X10  -  .6UiVJ  l 

EST2  1.06Q50  *  . 567X1 C  ,3 

Table  2  summarizes  the  relevant  statistic * 
for  these  specification  equations  the 


Figure  3.  One.Mlnutc  Average  Value  ef  ASVR1A 
XSVRIOO  and  RVR  lor  Advent  ion  Tec 
16-15  October  1978 
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Table  2.  STATISTICS  CATEGORY  I  SIKH  TAT!  ’■ 

Unlike  the  results  und<  r  ‘.aUg”*r\ 
11  conditions,  the  statistic*  of  Categ-:r>  : 
conditions  show  that  the  100  It  retrod  < EsT!  . 
out  performs  the  50  ft  method  (EST2/  b\  neari\ 
15  percentage  points  in  the  RMSE.  Similar 
differences  exist  m  the  other  statistics. 
Again  wt  find  the  remote  tower  cqu.iri,r*  j/e 
»upe<  mr  to  the  touchdown  RVR  m asurv-ent . 

To  demonstrate  these _results  w#  select 
ASVR200,  RVR  and  XTv*<2Qff  V redesignated  i >T1  i 
and  display  thtm  In  a  time  series  of  an 
advectlon  fog  episode  in  figure  6.  Note  the 
stubborn  bias  of  the  RVR  measurement  w  lie 
XSVR200  closely_jratt§  A5VR200. 
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Figure  4.  One  Minute  Average  Value-,  of 
ASVR200,  RVR  and  ASVR20&  for 
Advectlon  Fog  14-15  October  1978 

4.  PREDICTION  TECHNIQUES 

4.1  Markov  Model 

Gringorten  (1972)  adapted  for 
meteorological  use  a  special  class  of  the 
Markov  chain  called  the  Ornstein-Uhlenbeck 
process.  The  conditional  probability 
estimates  are  derived  from  the  uncategorized 
initial  condition  as  given  by  the  latest 
observation  of  the  predictand  and  are  based 
upon  an  exponential  decay  of  the  auto¬ 
correlation  coefficient  with  time.  An 
assumption  is  made  which  equates  the 
cumulative  frequency  distribution  (CFD)  of 
prevailing  visibility  to  the  CFD  of  SVR,  This 
Is  a  necessary  assumption  because  there  does 
not  exist  a  SVR  data  base  to  accurately 
determine  its  CFD*  A  twenty-three  year 
data  base  for  Otis  AFR  was  used  to  determine 
the  unconditional  CFD  of  prevailing  visibility 
as  a  function  of  time  of  day  and  season. 

SVR0  (SVR  at  time  tQ)  for 
Categories  1  and  11  conditions  were  obtained 
from  regression  equations  which  yield 
specifications  of  ASVR200  and  ASVR100.  These 
specification  values  were  used  as  inputs 
into  the  model* 

4 . 2  Regression  Estimation  of  Event 

Probabilities  (REEP) 

REEP  (Miller,  1964)  calculates 
probabilities  of  being  within  categories 
which  can  be  easily  converted  into  ex¬ 
ceedance  probabilities.  Instead  of  trans¬ 
forming  the  initial  SVR  into  a  most  probable 
SVR  from  which  an  exceedance  probability  can 
be  found  (Markov  model),  REEP  uses  the  Initial 
SVR  directly. 

Five  predictand  category  limits 
were  selected  to  coincide  wUh  the  thresholds 
prescribed  beforehand.  Category  limits 
were  assigned  to  each  predictor  based  on  the 
relative  frequency  distribution  of  the  SVR 
values  in  Che  dependent  data  set. 

REEP  is  formulated  to  insure 
internal  consistency  among  predictand 


The  classical  Markov  t  r  v*.s  1 1  io*> 
Matrix  M  can  be  used  In  preparing  forecasts 
for  any  number  of  time  steps  ini  :  .*  r* 

by  using  powers  oi  M.  A  prediction  (■■.*.; 
that  yields  probabilistic  forecast-  i-pr  r  « 
to  the  classical  Markov  process  but  w:i  .*  .* 
necessity  of  utilizing  M  explicits  \t, 

1968)  uses  the  coefficients  f  ron*  a  ■*. 

REEP  equations  to  determine  t  he  on-  -  i  •  y 
transition  matrix  P. 

Experience  with  the  model  resulted 
in  the  use  of  two  separate  transition  rat  rites 
to  cover  the  full  range  of  predictions.  For 
Category  It  operations,  a  five  minute  one- 
step  transition  matrix  computes  conditional 
probabilities  for  forecasts  ot  five  minutes 
or  greater.  For  Category  I  operations,  a 
two  minute  one-step  transition  matrix  is  used 
to  compute  two  and  five  minute  forecasts 
and  a  ten  minute  one-step  transision  matrix 
is  used  to  compute  forecasts  ot  ten  minutes 
or  greater. 

5.  RESULTS:  PREDICTION  TF  HN ! QlT ^ 

The  probability  forevu-l-  are 
evaluated  through  the  use  ol  the  P-m  >re 
as  defined  by  Epstein  1 1969)  >in(j  tfu rph\ l  l'>t 
and  by  reliability  graphs.  The  P-st  >r* 
represents  the  mean  squared  difference 
the  forecast  and  observed  probability  dis¬ 
tributions*  The  three  prediction  techniques 
were  applied  to  each  of  the  twenty-five 
independent  data  episodes.  The  results  of 
each  episode  were  combined  to  reflect 
the  overall  accuracy  and  reliability. 

Verification  results  of  the  f- 
score  for  Category  1  operations  given  in 
Table  3  reflect  that  the  Equivalent  Markov 
technique  yields  a  slight  improvement  over 
the  Markov  and  REEP  techniques  (2  p*  rcent 
and  9  percent  respectively).  Figures  to  7 
depict  the  reliability  oi  the  three  techni.;.-v- 
for  10  minute  predictions.  Clearly  tin* 
predicted  probabilities  oi  the  Equivalent 
Markov  and  REEP  techniques  are  qu i r <  close 
to  the  observed  frequencies.  However,  tie 
predicted  probabf  1  i  t  tes  of  the  Marl.ov 
technique  tend  to  underestimate  the  ob-vtwj 
f  requenctes. 
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Table  3.  P-SCORES  CATEGORY  1 


Table  A  suirntri/t  <•  *1  r< 
Category  11  operations*  Tin  i  ,  ■  v  ■  :*-• 
technique  yields  s\ stemat  i «  a  1 1 \  ‘  i  t  t  •  • 
over  the  Markov  and  KN.P  te-l.n:  : 
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Figure  3* 


Reliability  Graph  Markov 
Category  1 
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Figure  7.  Reliability  Graph 
Category  I 
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Table  A*  P-SCORES  CATEGORY  11 

A  dropot  i  in  skill  ol  I  he  Npiivale-t 
Markov  technique,  as  the  forecast  interval  N 
increased  (rom  10  to  bO  minutes,  is  due  in 
part  to  decreasing  resolution  in  the  low- 
probability  Corecasts*  The  percent aue  of 
forecasts  in  the  0-5  percent  range  lowers  iron 
85  to  72  to  58  present . 

6.  CONCLUSIONS 

Analysis  of  the  data  col  let  toil 
at  the  AFGL  WTF  demonstrates  the  accurate  anj 
reliability  of  a  remote  tower  system*  lor 
both  Categories  l  and  11  operations  such  a 
system  would  provide  a  probability  ot  de¬ 
tection  of  below-)  itni  t  conditions  ot  U1 
percent  or  better,  a  large  improvement  over 
using  the  surface  RVR  instrument.  The 
Equivalent  Markov  prediction  technique  vis 
shown  to  provide  accurate  and  re  Nab  V 
forecasts  ot  below-limit  SVR  conditions  and 
yield  slightly  better  results  on  indept ndrnf 
data  than  did  the  Markov  and  KEEP  techniques. 
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